Heavy quarkonium correlators at finite temperature: QCD sum rule approach 
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We investigate the properties of heavy quarkonia at finite temperature in detail using QCD sum 
rules. Extending previous analyses, we take into account a temperature dependent effective contin- 
uum threshold and derive constraints on the mass, the width, and the varying effective continuum 
threshold. We find that at least one of these quantities of a charmonium changes abruptly in the 
vicinity of the phase transition. We also calculate the ratio of the imaginary time correlator to 
its reconstructed one, G/Grcc, by constructing a model spectral function and compare it to the 
corresponding lattice QCD results. We demonstrate that the almost constant unity of G/Grcc can 
be obtained from the destructive interplay of the changes in each part of the spectral modification 
which are extracted from QCD sum rules. 
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I. INTRODUCTION 

In- medium property of heavy quarkonia provides infor- 
mation on the confinement-deconfinement transition in 
QCD. In relativistic heavy ion collisions, final yields mea- 
sured through dilepton channels depend on whether they 
can exist as bound states or not. If the deconfined plasma 
is produced, color Debye screening will melt the quarko- 
nia then suppress the resultant yields [Ij] . Although J ftp 
has been measured in heavy ion collisions at various ener- 
gies, no quantitative understanding has been reached yet, 
because of the intrinsic complexities of processes in the 
heavy ion collisions. See Refs. [2|-[3| for recent reviews. 
Therefore, it is important to investigate the properties 
of the quarkonia in an ideal environment to give a solid 
foundation, not only on the existence of a bound state 
but also on their detailed spectral modification such as 
the mass shift and broadening at finite temperature or 
density. In this respect, it was pointed out that a down- 
ward mass shift of J/^ in hadronic matter was caused 
by the decrease of string tension, which had been pre- 
dicted by lattice QCD, and as such can be a precursor 
phenomena of the confinement-deconfinement transition 
0. Therefore, the detailed determination of the spectral 
properties can play a key role in the study of QCD phase 
transition. 

Properties of the bound states have been traditionally 
investigated with quantum mechanical potential models. 
It is known that the mass spectrum of heavy quarkonium 
can be described well by the so-called Cornell poten- 
tial, which implements the Coulomb potential at short 
distance and linearly rising one at long distance @, 
This approach can be extended to finite temperature by 
assuming that all the effects of temperature can be ac- 
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counted for by the temperature dependent potential 0. 
In this approach, however, how to construct the potential 
relevant for Schrodinger equation is a non-trivial problem 
and various types have been examined by incorporating 
properties known from lattice QCD |9l-[l3l|. Although the 
potential model approach is related to QCD only through 
the temperature-dependent potential computed with lat- 
tice QCD, recent development in weak coupling methods 
such as pNRQCD |14l - [la | which is an effective field theory 
of QCD and a resummed perturbative approach p7l . [l8| 
shed a light on more rigorous foundation of the heavy 
quark potential at finite temperature. 

Direct evaluation of the quarkonium properties with 
lattice QCD has been carried out through the maximum 
entropy method (MEM) (l9l42l| . One can reconstruct 
the spectral function of a given channel by inverting the 
dispersion relation of the current correlation function cal- 
culated in the imaginary time. Indeed, Ref. indicated 
the existence of J/^ bound state even in the deconfined 
phase up to T ^ l.BTc. At high temperature, however, 
lattice QCD suffers from the limited size of the temporal 
direction, which also affects the accuracy of the recon- 
struction of the spectral function in MEM [2l|, . In 
this approach, the only reliable information seems to be 
the presence or disappearance of the first peak in the 
spectral density. A potential model calculation [l^ [23| 
indicates that the first peak observed in MEM at high 
temperature can be attributed to the threshold enhance- 
ment and that J/ip has already melted a.t T — 1.2Tc. 

Recently AdS/QCD approach has also been applied 
to the heavy quarkonium in medium [l^, HH] • Although 
a direct relation to real QCD is still missing, the ap- 
proach seems to give another insight to the problem from 
the viewpoint of the strongly coupled gauge field theory. 
Both Refs. [U, show notable spectral change in J/ip 
around and above Tc- 

In previous works [26l - l29t . we have proposed another 
approach to study the properties of heavy quarkonia at fi- 
nite temperature based on QCD sum rules; the a ppr oach 
extended the previous studies at nuclear medium j30l.f3l|. 
The QCD sum rule [H, H^l provides a systematic frame- 
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work which connects the current correlation function at 
deep Euclidean region to the spectral function integrated 
with respect to the energy variable with a weight that 
makes the integral dominated by the lowest pole. It 
has been applied to various aspects of hadrons quite 
successfully 0, [s^. Due to the asymptotic freedom, 
one can reliably compute the correlation function at the 
deep Euclidean region using perturbation theory via the 
operator product expansion (OPE) which provides non- 
perturbative correction through QCD condensates. For 
a heavy quarkonium, to a good approximation, one can 
truncate the expansion at the lowest dimensional local 
operator, which is the dimension four gluon condensate. 

The aim of this work is to extend our previous works 
to a more systematic analysis by incorporating the con- 
tinuum part of the model spectral function, applying a 
more sophisticated optimization procedure in determin- 
ing the spectral parameters, and then making a compari- 
son to the lattice QCD results. In this paper, we describe 
the detailed procedure based on the Borel transformation 
which is widely used in QCD sum rule applications. Then 
we discuss the spectral change of charmonia and bot- 
tomonia at finite temperature near and above T^. Using 
the spectral parameters obtained in the QCD sum rules, 
we construct model spectral functions and compute the 
imaginary time correlators, which we will compare with 
lattice QCD. 

This paper is organized as follows. In the next section, 
we briefly review the QCD sum rules for heavy quarko- 
nium at finite temperature, then explain the procedure 
based on the Borel transformation. We will show the re- 
sults of the spectral parameters in Sec. IIIII We will dis- 
cuss the imaginary time correlators reconstructed from 
the spectral parameters in Sec. IIVI Section |V] is devoted 
to the summary. 



II. QCD SUM RULES FOR HEAVY 
QUARKONIUM IN MEDIUM 

A. OPE for correlation function 

We start with the current correlation function 
n^{q) = i f d^xe"'-''{T[J^{x)J'\0)]). 



(1) 



We choose the currents for pseudoscalar (P), vector {V), 
scalar (S), and axial- vector (A) as 



j = ih-y^h 
f = hh 



(2) 
(3) 
(4) 
(5) 



with h being the heavy quark operator, c or h. For the 
axial-vector, we pick up the conserved components for 



Xc(b)i states. We will take the expectation value at fi- 
nite temperature. Therefore, in general, there are two 
independent components in both the V and the A chan- 
nels. We assume the momentum of the current to be 
q'' = ("^,0), i.e., a pair of quark and antiquark at rest 
with respect to the medium such that only one compo- 
nent becomes independent. Then we define the dimen- 
sionless correlation function as 



II^'^(<z2) 



n^'^(g) 



(6) 
(7) 



which can be expanded up to dimension four operators 
via OPE 



(8) 



with Go = {^G1^G<''^"') and G2 being the scalar and 
twist-2 gluon condensates, respectively. G2 is defined as 
the traceless and symmetric part of the gluon operator 



u^Ui, - ^gfj,^ ) G2 (9) 



which vanishes at T = according to the Lorentz invari- 
ance. The medium four velocity is set to (1,0,0,0). 
Hereafter we assume that all medium effects are imposed 
on the change of the local operators [s^ . This assump- 
tion is justified when the typical scale of the condensates 
is smaller than the separation scale [2^, namely. 



Ami - 9' > (G) - (Aqcd + aT + by.f 



(10) 



with TO/i = rric or mb being the heavy quark mass. The 
large heavy quark mass mh allows us to work even at 



the physical energy scale q'^ = nij , to|- and so on, al- 
though somewhat marginal in reality [s^ . In this case, 
the OPE gives a formula for the bound state mass that 
is proportional to the change of the color electric field 
squared. This formula is the QCD second order Stark 
effect [srl-lsg^ that can also be obtained from the leading 
order correction of the static potential in pNRQCD pi| . 
Combining the formula with the temperature dependence 
of the electric condensate, one finds a sudden mass shift 
at Tc [39]. In the QCD sum rule, we go to the deep 
Euclidean region = lj^ = —Q^ ^ 0, in which the 



condition pOj) is well fulfilled. Therefore the Wilson co- 
efficients Gi are the same as the vacuum case and have 
been calculated in Refs. [H, HO, IH, . See Ref. [H for 
a list. 
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B. Borel transformation and dispersion relation ,— , 



The Borel transformation of the correlation function is 
defined as 



iQ 



2\n+l, 



d 



M'\M^) = lim 

(11) 

If one does not take the limit, the derivative of the corre- 
lation function corresponds to the moment of the corre- 
lation function which was used in the previous sum rule 
works for the heavy quarkonia [26l - [28l Isol . [sil Issl . |40| . 
Taking this limit corresponds to going to deeper Eu- 
clidean region for better perturbative expansion while 
retaining the connection to the resonance through large 
n [l^l- Indeed, Eq. pUj) is expected to be better satis- 
fied as the typical OPE term (^^^^_„2\d (C*^) reduces to 



-4ml /M^ 



(G"*) after the transformation, and hence 
the condensate contribution is further suppressed by ^. 
In the heavy quarkonia, the moment sum rule works well 
enough to extract the mass due to the large separation 
scale coming from the heavy quark mass. Nevertheless, 
the Borel transformation approach has several advan- 
tages for more systematic analysis as revealed below. 

For the expanded heavy quarkonium correlation func- 
tion the Borel transformation can be analytically car- 
ried out as 



+ c\v)<j>,{T)] 



(12) 



with a dimensionless scale parameter v = Am\/M'^ . The 
first line of Eq. is the same as that derived in 

Ref. [4l| except for the temperature dependency of the 
scalar gluon condensate term cjjb- The second line shows 
the twist-2 term which appears in the case of medium. 
4>b and (pc are defined as 



47r2 



9(4m2)2 

47r2 
3(4^ 



Go(r), 



G2(r), 



(13) 
(14) 



as given in Ref. [2l[23l- While A\v), a'^ {v) and ^ {v) 
are given in Ref. [41|, the transformed twist-2 coefficient 
{v) is derived for the first time in this paper. For com- 
pleteness, we list all the Borel transformed Wilson coef- 
ficients used in Eq. ([T2|) in Appendix \^ While Bertl- 
mann worked on the on-shell renormalization of heavy 
quark mass in Ref. j41|], we maintain the off-shell renor- 
malization as a straightforward extension from Ref. [i^ . 
Hence, the correction term —^^^h'^iv) in a'^ {v) is in- 
cluded throughout this calculation. Note that this term 
is also necessary to keep the perturbative correction term 

small enough in A^(Af^) and in —^^fjjjij- that is used 
later. Since this part is temperature independent, the 
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FIG. 1: (Color online). Temperature dependence of the gluon 
condensates Go (T) and G2 (T) (top) and effective temperature 
dependent coupling constant (bottom). 



difference does not affect our aim but enables us to pro- 
ceed in a more transparent way by retaining the rela- 
tion with the previous moment sum rule analyses. In 
Eq. ()12|) , external inputs are heavy quark mass mh (con- 
tained in A{v), ipb and 0c), strong coupling constant 
as{M^) and gluon condensates Go{T) and G2[T). In 
this paper, we put mc{p^ = —m^) = 1.262 GeV in- 
ferred from the rndp^ = —2m^) = 1.24 GeV used in 
the previous works [26|-|30| and mb{p^ = —rnl) ~ 4.12 
GeV. as(Af^) is calculated from the running coupling 
formula from as{^m'^) = 0.21 also used in the pre- 
vious works. as{9>rnl) = 0.158 is used for bottomo- 
nia. The gluon condensates have been extracted from 
the results of pure SU(3) lattice gauge theory [i^l and 
shown in Fig. [T] Here effective temperature dependent 
coupling constant as{T) is used for the determination 
of G2(T) which is the symmetric and traceless part of 
the gluon operator. This is done with the identifica- 
tion i^STiG'^^^G;''')) ^ ^{STiGl^G'}^^)) based on 
the the separation scale in the OPE. We adopted aqq{T) 
shown Ref. [isl and took the value at r = rscrccn. It is 
also shown in the bottom panel of Fig. [TJ 

The current correlation function is related to the spec- 
tral function through the dispersion relation. At finite 
temperature, the spectral function is given by the imag- 
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inary part of the retarded correlation function 11^ (q) 
which is in general different from Eq. ([Ij. However, we 
can relate it to the spectral function by virtue of the fact 
that n^(cj) = n(a;2) since oj^ = -Q^ < [se^. Putting 
tanh(y^/2T) = 1, which is safely satisfied for ^ TO/j/i, 
we have the dispersion relation for the Borel sum rule in 
the same form as the vacuum case 



Jo 



(15) 



Note that the weight factor in the dispersion integration 
is now exponential e~*/*^ while it was the inverse power 
(s + Q^)^" in the moment sum rule. 



C. Analysis procedure 

Assuming the quark-hadron duality, we take a model 
spectral function of a given current as a simple ansatz for 
the imaginary part of the correlation function and call it 
the phenomenological side. First we decompose it into 
the pole and the continuum contribution as 



Imn(s) = ImnP°"=(s) + Imn™"*(s), 



(16) 



with 



ImnP°'<=(s) 
Imff°"'(s) 



foSis - m^), r = 

, r > 0, s> 4ml 



(s — TO^)^ + sF^ 

(s-so)ImII'^'P"''(s) 



(17) 
(18) 



The pole term is the same as in the previous works ^2^- 
[28| . We consider possible finite width in the deconfined 
medium, because decay into hh pair, which was forbid- 
den in vacuum due to the Okubo-Zweig-Iizuka (OZl) rule, 
becomes possible. This is done by implementing a rela- 
tivistic Breit-Wigner function and cutting off the contri- 
bution below hh threshold to avoid possible numerical 
artifacts as discussed later in Sec. IIIIBl 

We adopt the perturbative part of the spectral func- 
tion including as correction but with the sharp threshold 
factor 9{s — sq) as a model for the continuum; such form 
reproduces the corresponding part of the OPE side when 
putting So = 4m^. These functional forms are explicitly 
given in Ref. 3j| and listed in Appendix [B] for the com- 



pleteness. Since there are known excited states such as 
-0' between the lowest lying state and the physical con- 
tinuum threshold, one may think this model is an over- 
simplification of the real spectrum. However, due to the 
suppression coming from the Borel transformation, this 
simplification does not affect the property of the low- 
est pole. Instead, this form results in a little smaller 
continuum threshold value than that from the analysis 
incorporating the excited states explicitly as we will see 
later. Moving the continuum part to the OPE side in 



Eq. (|T5|) . one can isolate the pole term. There is an ad- 
ditional contribution to the spectral function from the 
absorption of the current by the thermally excited par- 
ticle, i.e., Landau damping which shows up as a peak at 
s — when qr = 0. This was recognized in Ref. [4J] in a 
QCD sum rule framework and has been called the "scat- 
tering term" . Recently it has been emphasized that this 
gives constant contribution to the imaginary time corre- 
lator which is the basis of the spectral function study in 
the lattice QCD [4^. In the QCD sum rule application 
in the deconfined phase, we can neglect this contribu- 
tion, as explained in Ref. [1^ [2^ . The scattering term 
appears in the OPE through the bare heavy quark con- 
densates (hVDD..Dh) which is converted into the gluon 
condensates via heavy quark expansion at T = [46j . 
To a first approximation that assumes free heavy quarks 
in a medium, we can put the same quantity on the phe- 
nomenological side so that it cancels the corresponding 
contribution in the OPE. 



Differentiating both side of Eq. ([15]) with respect to 
1/A/^ and taking its ratio to the original equation, one 
has 



d 



5(1/M2) 



[A1(m2) -X'=°"'(m2)] 



I dsse-"/^^'lmnP°'°(s) 

Jiml 

h 

/•OO 

/ dse-"/*^'lmIIP°^'=(s) 



(19) 



where A^'^°"'(M2) is the Borel-transformed continuum 
spectral function according to Eq. (jl5l) . One immedi- 
ately finds that the right-hand side of Ec^. ([T9l) gives the 
squared pole mass for F = 0. There are three spec- 
tral parameters to be determined in Eq. pole mass 
m, width F and effective continuum threshold sq. The 
strength parameter /o or / contained in the pole term 
cancels by taking the ratio. We solve Eq. (IT^ for the 
mass m as a function of the Borel mass M"^ (which we 
call Borel curve) with given sets of F and sq- While the 
extracted mass m depends on through the left-hand 
side of Eq. by construction, it should not do so since 
is an external parameter. The apparent dependence 
of m on is due to the truncation of the OPE and the 
insufficient subtraction of excited states and the contin- 
uum part in the spectral density. The truncation of the 
OPE shows up as strong dependency at small 
while the insufficient subtraction does at large M^. In 
practice, however, we can expect independent m at 
intermediate region after tuning the effective thresh- 
old parameter sq in the continuum part [s^ |43| . Intro- 
ducing the finite width also affects the m{M'^) in the 
small A/2 region [4^. We thus introduce the following 
quantity to determine the best set of the parameters that 
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gives the flattest curve in the intermediate region; 



,2 _ 



1 



M2 

max 



Mi 



Mi 



2\i2 



dA'rimiA/r) -m{M^)] 



where is defined by dm{AP)/dM'^ 



(20) 
The 



:A^max] is called 



Q in uciiiicu uy uuiyin j / uiu \m^=AI^ — 

range of the intermediate S [Af^^ 
Borel window, in which the convergence of the OPE and 
the pole dominance of the dispersion integral are satis- 
fied. We fix M^jjj by requiring the dimension four cor- 
rection terms to be smaller than 30% of the total since 
it is expected that the contribution from the next higher 
dimensional operator is kept less than 10% of the total 
within this condition j32j . For bottomonium systems, 
while this condition is always fuUfiUed in typical val- 
ues due to the larger quark mass in Eqs. and the 
perturbative radiative correction term as{M'^)a'^ {f) can 
be large enough to spoil the perturbative expansion. We 
therefore impose this term to be less 0.3 at > M^j^. 
The pole dominance should be also imposed on the cri- 
teria to preserve the reliability on the extracted property 
of the pole part of the spectral density. We determine 
M^g^x by requiring the continuum contribution to be less 
than 30% of the total perturbative term. As we shall see, 
reducing the value of the continuum threshold makes the 
Borel curve flatter. Therefore, this criterion does not af- 
fect the pole mass at the best fit while it becomes impor- 
tant when the Borel window is narrow. The measures 
the average deviation of m(M^) from its value at the A'P 
"plateau" characterized by diTi{M'^) / dM'^\]^j2^]^,j2 = 0. 
While it vanishes in the case of completely indepen- 
dent mass, can be regarded as a systematic error on 
the extracted mass ui^Mq) such that m(MQ ) -I- a/x^ gives 
the upper limit. Note, however, that it gives much larger 
deviation than actual uncertainty when Borel window 
contains the strongly dependent part of the Borel 
curve at small which is due to the truncation of 
OPE. Hence, a refinement on the determination of the 
Borel window will give more quantitative insights on the 
uncertainty in the extraction of the mass parameter. 

To illustrate the minimization procedure, we show be- 
low how the Borel curve obtained by solving Eq. 
changes with respect to the external parameters. We de- 
pict some examples in Fig. [5] Let us start with the dot- 
ted curve corresponding to T = 0, ^/sq = 3.5 GeV, and 
F = for J/ip- ^From the definition of the Borel mass 
JVP , one sees that the Borel curve looks like the n depen- 
dence of the moment ratio shown in Refs. [1^ [l^ except 
for the reversed direction of the horizontal axis. Here 
we draw the lines only within the Borel window. There- 
fore the dotted line is truncated at AP = 1.12 GeV^ and 
3.69 GeV^. The arrow on each line indicates the loca- 
tion of Mg and the upper limit evaluated from \f~)^ ■ As 
temperature increases, the gluon condensates decrease as 
shown in Fig. [TJ This is reflected by the lowered curves 
obtained for T = 1.04Tc, as in line with the moment 
sum rule case [1^ [13]. Since the condensate contribu- 
tion becomes dominant in the OPE side at small A'P, 
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> 
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FIG. 2: (color online). Some examples of the Borel curve 
taken from J/^p at T = 1.04rc. Solid line shows the curve 
for y/so — 3.5 GeV and F = 0. Long-dashed line shows the 
case in which is reduced while in short-dashed line F 
is increased to 100 MeV. For reference, the thin dotted line 
shows the curve at T = with y/so — 3.5 GeV and F = 0. 
The arrows accompanied by a short horizontal line indicates 



the upper limit of the mass at M — Mq evaluated by 



3.2 



3.1 



> 
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7/i|/,T=1.12Tc 



So =3.5 GeV, r=0 
sj'^=3.08 GeV, r=274 MeV 
So ^=3.5 GeV, r=336 MeV 
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FIG. 3: (color online). Same as Fig. [1 but for T = 1.127^ 
See text for description. 



one sees great reduction of mass in this region. Then 
the solid curve, for T — 1.04rc, shows a minimum at 
smaller Borel mass. This corresponds to the shift of the 
minimum of the moment ratio to large n in the moment 
sum rule. Regarding the solid curve as the base line, one 
sees that decreasing the continuum threshold flattens the 
curve at large AP. In this case, however, reduction of the 
continuum threshold makes the AI^^^ smaller and thus 
the resultant Borel window becomes narrower. This is 
why the dash line ends at M^^^ = 2.08 GeV^. One 
also sees that the curve is almost flat above = 1.5 
GeV^ . This means that further reduction of the contin- 
uum threshold breaks the stability of the Borel curve, 
i.e., the mass decreases monotonically and no minimum 



6 



would exist. M^^^ does not change against the reduc- 
tion of the continuum threshold since it depends only on 
temperature through the power correction terms in the 
OPE. On the other hand, if one increases the width, it 
raises the mass, especially at low as clearly seen in 
the short-dashed curve. One sees that, at AP far from 
the Af^jjj, the two lines, one obtained by reducing the 
continuum threshold and the other by introducing the 
width, show almost similar flatness. In the present case, 
the rapid rise in the T — 100 MeV gives ~ 8.55 x 10"^ 
GeV^ which is much bigger than = 1.89 x 10~^ GeV^ 
of the ^/so = 3.2 GeV curve. However, this depends 
on the choice of the criterion in the determination of the 
Borel window. If one tightens the criterion, to 10% power 
correction for instance, M^^^ becomes larger and then 
of the r = 100 MeV will be smaller. This indicates the 
difficulty in accurately determining the spectral parame- 
ters when one takes into account the change of both the 
width and the continuum threshold. Nevertheless, one 
can make the curve flatter by decreasing the continuum 
threshold without introducing broadening up to a certain 
temperature. Note that the arrows become shorter as the 
curve becomes flatter. As explained above, while there is 
about 50 MeV uncertainty in the largest case, it is due 
to the strong dependence seen in small region. 
Therefore one should not take these values so seriously. 

This situation changes if one goes to higher temper- 
atures. We plot some examples from T = l.l2Tc in 
Fig. [3l As the solid line shows, no stability is achieved 
in ^/so = 3.5 GeV and F = case. This is similar to 
what is seen in the moment sum rule above T > 1.05Tc 
[26j . Now we can try to restore the stability by varying 
So and F. From what we learned from Fig. [51 reducing 
So decrease m(M^) especially at high M^. In this case, 
however, Borel window closes before stability is restored; 
-^^max < -^min occurs. Therefore, one has to increase 
width to recover the stability. In other words, the break- 
down of the stability occurring above Tc can now be re- 
garded as the onset of the broadening. We denote this 
temperature as Tonsct, which depends on the channel as 
we shall see below. Note that this does not mean F must 
be below Tonsot, since one can find the best parameter 
set with F > at T < Tonsct after an additional con- 
straint is given. Tonsot should be regarded as the upper 
limit of temperature at which broadening sets in. One 
should also note that Tonsct depends on the criteria for 
the Borel window. One can broaden the Borel window 
by relaxing either or both of the criteria. For example, 
if one sets the continuum contribution to be less than 
50% instead of 30%, A^^ax becomes larger thus can open 
the Borel window. Indeed such a situation is realized in 
some cases considered in this paper, as seen in the resul- 
tant spectral parameters summarized in Tables IVHIXIIll 
When the x^ takes its minimum at the smallest ^/sq sat- 
isfying the criteria, making the Borel window larger by 
relaxing the M^^x criterion leads to smaller ^/sq while 
retaining the width. Note, however, that there might be 
"pseudopeak" artifact in the Borel curve for a too relaxed 



criterion [49| . 

The long-dashed line in Fig. |3] denotes the case in which 
we introduce F — 274 MeV with decreasing the contin- 
uum threshold to ^/sq ~ 3.08 GeV. As a result of raising 
the Borel curve at small M"^ while lowering it at large 
Af^, the shape of the curve becomes convex contrary to 
the lower temperature cases. As seen in the short-dashed 
line, one can restore the stability only if one increases the 
width. If we keep ^/s^ = 3.5 GcV, the resultant width 
is 336 MeV. The values of x^ are 1.69 x lO''^ GeV^ for 
y/s^ = 3.08 GeV and F = 274 MeV and 1.9 x 10"^ GeV'^ 
for for y/so = 3.5 GeV and F = 336 MeV, respectively, 
indicating the almost equally flat curves and again the 
difficulty of comparing the curves by varying both sq and 
F. In Fig. 131 a/x^ for the two stable Borel curves are a 
few MeV, which are small enough to be neglected in the 
figure. 

In the following, we use the evaluation using 
Eq. (j20|) only for determining the best curve among those 
cases with the same continuum threshold but with dif- 
ferent F, in order to avoid the biases imposed by the 
choice of the Borel window. In this way one fixes one 
edge of the curve M^^x so that x^ measures only the 
effect of introducing the width. For some cases where 
F = always gives the flattest curve, we may use x^ to 
determine the best sq value, since the other edge of the 
Borel curve is fixed. For example, we can safely deter- 
mine the best sq by evaluating x^ at T = 0. Furthermore, 
though we maintain the criterion for the Borel window 
as explained, we can easily estimate how the best value 
changes with respect to the change of the criterion. If 
one relaxes the pole dominance condition, it extends the 
Borel window to larger Af ^ therefore continuum thresh- 
old giving the best x^ will become smaller. On the other 
hand, if one requires the smaller power correction, M'^^^ 
becomes larger and thus the x^ will be more sensitive to 
the continuum. As long as we preserve reasonable values 
of these criterion, typically 10 — 30% for power correc- 
tion and less than 50% for the continuum contribution, 
we find the uncertainty of the obtained mass to be about 
a few tens MeV. Since we maintain the same criterion 
even at different temperatures, the relative in- medium 
change of the spectral parameters is not affected by the 
particular choice of the criterion for the Borel window. 

Before closing the section, we would like to address 
possible uncertainty on the extracted parameters. The 
source of the uncertainty is roughly classified into two 
parts; one is temperature independent and the other is 
dependent quantities. The former consists of the giuon 
condensate at T = 0, strong coupling constant, and the 
heavy quark masses. Whereas these quantities certainly 
affect the value of the spectral parameters, the relative 
changes at finite temperatures from the vacuum values 
do not differ by changing them within the constraints 
from the experiment. Hence we focus on the effect of the 
temperature dependent part here. We estimate the la 
and 2a uncertainty of the OPE side through the temper- 
ature dependent part of the gluon condensates by reading 
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5m,/^ [MeV] 

FIG. 4: (color online). Constraint on the mass shift and 
the width of J/V" at T = 1.04rc and ^ = 3.28 GeV. In 
addition to the thick lines denoting the constraint, contours 
of the equal values (unit of deviation from the minimum) 
are plotted as thin lines. 




FIG. 5: (color online). Same as Fig.|4]but for = 3.5 GeV. 



off that of the energy density and pressure in the lattice 
results shown in Ref. [42|. Then, we extract the spec- 
tral parameters with the OPE side shifted by ilcr and 
±2(7. The resultant constraints which show up the typi- 
cal correlation between the mass and the width together 
with the uncertainty deduced from the deviation of the 
OPE side are displayed in Figs. 2] and [5] The solid line in 
Fig. |4] shows the case corresponding to the long-dashed 
line in Fig. [2] the best Borel curve is obtained by reducing 
the threshold parameter without introducing width. This 
fact is reflected to the kink of the Xmin contour (thin solid 
line). If the OPE side increases, the minimum shifts to fi- 
nite r (see next section). The obtained constraint shows 
a clear correlation between the mass and the width as 
pointed out by us in Ref. 26]. Quantitatively the rela- 
tion can vary within ±(10 — 20) MeV due to the change 
of the temperature dependency. Nevertheless the overall 
behavior of the correlation is preserved. Figure [S] shows a 
similar case but the finite width ~ 100 MeV was chosen 



0.3 




1 1.5 2 2.5 3 3.5 4 
Borel Mass [GeV^] 



FIG. 6: (color online). OPE coefficients of V (upper panel) 
and S (lower panel) channels in charmonium systems. Solid 
lines denote the temperature independent perturbative cor- 
rection term asiM^^aiv). Long and short dashed lines stand 
for the scalar gluon condensate term b{i')<f)f,{T) for T = and 
T — Tc, respectively. Dotted line shows the twist-2 contribu- 
tion c(j/)0c(r)at T = Tc. 



as the best fit. The correlation between the mass shift 
and the width does not differ from the other example 
but the shape of the contour slightly does, refiecting 
the straight contour of Xmin- these figures, the area 
surrounded by the la or 2a correlation lines and one of 
the contour can be regarded as a possible region of 
the mass shift and the width including uncertainty for a 
given 



III. TEMPERATURE DEPENDENCE OF 
SPECTRAL PARAMETERS 

A. Implication from OPE side 

Let us begin with examining how temperature depen- 
dence of the gluon condensates is related to the spectral 
parameters by looking at the dispersion relation (jlSp . 
The temperature dependence of the OPE side [Eq. p^ ] 
comes from the gluon condensate terms b4>b and c0c- In 
Fig.ini we plot the OPE coefficients in Eq. ((T^). One sees 
that both the scalar gluon condensate contribution and 
the twist-2 one increase as temperature increases. These 
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dependencies come from the fact that Go{T) and G2(T) 
are monotonically decreasing functions of the tempera- 
ture while the coefficients b{i') and c(i/) have always neg- 
ative sign.^ Therefore, up to dimension four operators, 
M{AP) in the OPE side is always increasing function of 
the temperature. 

One also sees that the expansion coefScients of the P- 
wave state are larger than ^-wave's at the same Borel 
mass M^, as also seen in the moment sum rule case 
[H, [2^ ■ These properties result in the correct mass split- 
ting between the J/ip and the Xc states, and induce larger 
mass shift for the P-wave states through the derivatives 
of these coefHcient with respect to From the be- 

havior of the condensate contributions at low AP, one 
realizes how the location of the Borel window in the case 
of charmonium changes with respect to the temperature. 
At temperatures close to Tc, M^-^^ is determined by the 
scalar condensate and thus becomes smaller as temper- 
ature increases. It eventually starts to increase when 
the twist-2 contribution dominates over the scalar one, 
|c(z.)0,| > \bi,y)q^b\. 

^From the dispersion relation (fTS]). where the weight 
of the integral over the spectral function is positive defi- 
nite, one obtains the constraint equation for the changes 
of the spectral parameters against the change of the OPE 
side discussed above. The phenomenological side (|16p- 
([T5)) has four parameters; effective continuum threshold 
So, pole mass m, width F, and overlap /. If only one of 
these four quantities is allowed to vary as temperature in- 
creases, the respective changes of the parameters needed 
to match the OPE side are, 

• So ■ decrease, 

• m : decrease, 

• F : increase, 

• / : increase. 

In practice, all of these quantities can change not only 
to the expected direction but also to the opposite, as 
long as the total combined change of the spectral func- 
tion matches the OPE side. The Borel transformation 
procedure explained in the previous section provides an 
optimization way to find out the best set of the changes. 



TABLE I: Spectral parameters of cc systems at T = 0. Ex- 
perimental masses are taken from Particle data book [59| . 

System ^i^ [GeV] m [GeV] mexp [GeV] [GeV^] /p [GeV^ 



J/i' 

Xcl 



3.48 
3.54 
3.82 
3.78 



2.993 
3.060 
3.406 
3.470 



2.980 
3.097 
3.415 
3.511 



1.547 
1.971 
2.552 
2.810 



0.396 
0.393 
0.303 
0.196 



3.20 
3.15 
3.10 
3.05 
8 3.00 
^ 2.95 
2.90 
2.85 
2.80 
0.35 
0.30 
„ 0.25 
I 0.20 
" 0.15 
0.10 
0.05 
0.00 



T=0.87Tc 
0.98T,. 
I.OOTj, 
1.02T^ 



-X- 



1.047^ 
1.077^ 
1.097^ 



-O- 



7* 



3.1 3.2 3.3 3.4 3.5 3.6 



[GeV] 



FIG. 7: (color online). Constraint on the spectral parameters 
obtained from the QCD sum rule for rjc- Upper and lower pan- 
els show the extracted masses and widths as functions of the 
continuum threshold, respectively. Symbols guided by lines 
represent the different temperatures. The vertical and hori- 
zontal lines indicate the continuum threshold and the mass at 
T = 0, respectively. 



B. Results for charmonia 

We carried out the analyses for the charmonia for var- 
ious temperatures around Tc- For a reference, we sum- 
marize the results of T = in Table HI One sees that 



^ In P channel b(u) can be positive, as understood from Eq. ||A2TJ, 
but we found it is mostly negative for values of u corresponding 
to Mq and M{M^) retains the property of increasing function 
of the temperature. 



the masses are well reproduced by the common parame- 
ter set indicated before. Finer tuning on the quark mass, 
the coupling constant, and the gluon condensate may im- 
prove the small discrepancies with the experimental data 
but we are not intending to do so here since our aim is 
to investigate the relative change from the vacuum value 
induced by the medium. 

We display the charmonium masses and widths ex- 
tracted using the method described in the subsection III CI 
as functions of the effective continuum threshold ^/sq in 
FigslTlfTOl One sees the constraints among m and F 
for various temperatures around Tc- Qualitatively all the 
four cases are quite similar; the mass increases almost lin- 
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FIG. 8: (color online). Same as Fig. [T] but for J/ip- 
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FIG. 9: (color online). Same as Fig. [T] but for XcO- 



early with increasing ^/sq at all the temperatures. While 
it shows linear increase, the slope of the mass changes 
depends on whether the width is zero or not. One can 
see kinks in the mass curves at the same values of the 
horizontal axis as those of the width, which vanishes at 
small Although we do not adopt a/x^ as a way to 

determine the best ^/sq, we notice that it takes the min- 
imum near the kink among various ^/so values at a fixed 
temperature. At T = TLnsotj beyond which width must 
be nonzero to maintain the Borel stability, the widths 
shown in the lower panels show linearly increasing be- 
havior with increasing continuum threshold. The onset 
temperatures are summarized in Table |lTl One sees rjc 
starts to broaden earlier than J/ tjj while Tonsot of the P- 
wave states are the same. This may indicate a different 
temperature effect on the spin-spin interaction responsi- 
ble for the mass splitting in the iS-wave states. In com- 
parison with what we learned in Sec. IIII Al one sees the 
result is more complicated than the simple analysis by 
the dispersion relation. One sees the reduction of the 
mass always couples with that of the effective threshold. 
This is a consequence of the optimization by the Borel 
transformation, i.e., using Eq. (jl9p and looking for the 
stable curve. To see the reason more explicitly, one can 
go back to Fig. [2] The reduction of the condensates low- 
ers the mass without any change in other parameters (see 
dotted and solid curves). The requirement of the Borel 
stability makes the curve flatter by reducing the effective 



threshold (see dashed curve). Thus the downward mass 
shift always occurs with the reduction of the effective 
continuum threshold if there is no broadening. 

Behavior of the width at small ^/sq and at T > Tonsct 
seems different from the genuine linear behavior. It first 
decreases as ^/sq increases, then turns to increase. For- 
mally we could obtain stable Borel curves at higher tem- 
peratures than those shown in the figure. At such high 
temperature, the width is always at an order of hundreds 
MeV and the Borel curve is similar to those displayed in 
Fig. [3] We would like to stress, however, that this re- 
sult may not be a physical one; in this region, the mass 
is also small while the width becomes 100 — 200 MeV 
or more. Clearly the Breit-Wigner ansatz in the phe- 
nomenological side ([T71) which cut off the lower energy 
part than 4m^ does not match with the dispersion in- 
tegral dTSl) . If we do not impose the cutoff, the strong 
suppression factor combined with the Breit-Wigner form 
in the Borel-transformed dispersion relation (|15l) leads 
to numerical artifacts such that the contribution com- 
ing from spectral function much below 4m^ comprises 
a subdominant fraction of the total dispersion integral. 
We depict an example taken from J/tp a.t T = l.OTTc 
in Fig. [TTJ The same consideration also holds for ijc, 
XcO, and xa at T > l.OTTc- One sees that the inte- 
gral receives large contribution from the energy region 
much smaller than 4m^ when the width becomes larger, 
despite the increase of the mass. For example, let us 
consider the change of the width when lowering the hh 
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TABLE II: Onset temperatures of the width Tonsct for char- 
monia 



3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 
[GeV] 

FIG. 10: (color online). Same as Fig. [T] but for Xci- 
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FIG. 11: (color online). Integrand of the dispersion integral 
(|15|l with the pole term obtained in T = l.OTTc. We normal- 
ized the different lines so that they become unity at s = rr? . 
The thin dotted line parallel to the vertical axis indicates the 



threshold by 1 GeV^ in the data shown in Fig. [TTJ With 
this change, F = 64 MeV at ^ = 3.2 GeV and F = 110 
MeV at ^ = 3.4 GeV become 30 MeV and 62 MeV, 
respectively. This is so because the contribution to the 
dispersion integral from s = 4m^ — 1 to s = 4m^ is large 
enough to compensate the smaller Brcit-Wigner width. 
We notice, however, that the solution of the sum rule, 
Eq. ([TO)) , does not exist at near M^j^^ for even smaller 
threshold as in the case shown in Ref. [Zs^. 

We also notice that this artifact is absent in the mo- 
ment sum rule up to n = 20 beyond which it breaks 
down. Therefore, changing lower limit of the integration 
range from 4rn^ to will not affect the previous results. 
At present, use of the vacuum dispersion relation, which 
cuts off the contribution below 4to^, seems effective to 
estimate the width when its magnitude is less than 100 
MeV. To give more quantitative results, we may need 
to take into account more detailed structure beyond the 
Breit- Wiener ansatz. Recent resummed perturbative cal- 
culation [50'| might provide useful information for a bet- 
ter modeling. Furthermore, Borel curves at low A/^ will 
be more influenced by higher dimensional operators we 
have neglected. Since the width is sensitive to the low 
behavior of the Borel curves, as shown in subsection 
III C[ it may receive sizable correction from those opera- 
tors. At present, temperature dependence of the higher 
dimensional operators is poorly known. More quantita- 
tive analysis of the width in the non-perturbative manner 
thus needs further efforts. 



Since Figs. ITHTOI give only constraints, one needs to 
specify one of those spectral parameters to discuss spe- 
cific temperature dependencies of each parameter. Pre- 
vious analyses |26l - [29l | correspond to sq — >■ oo hmit. For 
instance, if ^fs^ retains the vacuum value, the mass de- 
creases rapidly until T = Tc in r]c and T = 1.02rc in J/V"- 
We do not show the results of constant for the P- 
wave states, as Eq. (fT9|) has no solution at certain region 
inside the Borel window and thus [Eq. (UHl)] cannot 
be evaluated before it reaches the minimum as a func- 
tion of the width. This absence of the solution actually 
occurs in S'-wave cases also, especially at larger y'so and 
comes from the non-monotonic behavior of Borel trans- 
formation of the Breit-Wigner function 48|. Note that 



this does not mean the corresponding parameter sets are 
completely excluded, since one may choose another (nar- 
rower in most cases) Borel window such that the solution 
exists in any G [-^min; -^max]- Nevertheless, from the 
almost linear dependence of the mass and width on 
one can extrapolate the lines up to the desired value to 
have a rough estimate. Then, one finds in all the chan- 
nels that the mass first decreases then the width starts 
to increase, as the temperature increases when ^Jsq is 
held fixed. Note that this transition of the temperature 
dependence of the mass is caused by the start of the 
broadening of its width. We would like to point out that 
the analysis at sq = constant is not the same as that at 
So — ^ oo in which the determination of the flattest mass 
curve by does not make sense. When the mass is held 
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T/T, T/T, 

FIG. 12: (color online). Temperature dependence of the spec- 
tral parameters of J/ip (closed symbol) and 77c (open symbol) 
extracted from QCD sum rules combined with the second or- 
der Stark effect. The masses (top-left panel) are obtained 
from the second order Stark effect and used as inputs for the 
sum rule analyses. Other quantities are obtained by reading 
off the result of Figs. [7] and [8] which match with the mass 
shifts. 



to its vacuum value, the constraint is satisfied by the in- 
crease in both and the width. For example, in the 
J /ip case, ^/sq becomes 3.6 GeV and F = 128 MeV at 
T = 1.04rc while F = is still possible if the mass and 
decrease to 2.92 GeV and 3.16 GeV, respectively. 

At present, the temperature dependence of the con- 
tinuum threshold is not clearly known yet. In fact, our 
threshold parameter should be regarded as an effective 
one since we do not take radial excited states such as 
Tp' {2S) and r]c{2S) into account nor temperature depen- 
dent behaviors near the threshold [12,, J^, ^23|| . Neverthe- 
less, if those states dissociate at as low temperatures as 
Tc, one may regard as a physical threshold within 
temperatures between Tc and a certain temperature at 
which the pole position becomes so close to the threshold 
that the pole and the continuum part of the model spec- 
tral function starts to overlap. For the P-wave states, the 
model will be better due to the absence of excited states 
below the threshold. One might be able to interpret the 
asymptotic value of the quark-antiquark potential as the 
continuum threshold [11], which decreases as tempera- 
ture increases irrespective of the choice of the potential 
0, mi ■ This fact might be related to the decrease of the 
mass and subsequent dissolution of the D mesons [s^. 
In this case, the obtained constraints give the down- 
ward mass shifts in the all channels. If the reduction 
is strong, only the mass shift occurs without broadening 
up to T = Tonsct- If not, the widths will start to broaden 
gradually together with the moderate downward shift of 
the mass. 



Another external constraint can be obtained from the 
second order Stark effect in QCD [IMIIEI. Although 
the applicability to the charmonium systems is marginal, 
it gives a genuine downward mass shift due to the rapid 
increase of the color electric condensate |39i] . For illustra- 
tion in the case of downward mass shift, we combined the 
result of Ref. [s^ with those of Figs. [7] and HI by finding 
the masses in Figs. [7] and [S] that matches with the re- 
sults of the second order Stark effect and then looking at 
the corresponding continuum threshold and width. The 
results of the masses, continuum thresholds, binding en- 
ergies defined by i?bin = ^/sq — m, and the widths of 
J/ip and rjc are displayed in Fig. [121 As explained, the 
results for T > Tonsct are marginal. Moreover, the sec- 
ond order Stark effect has also limitation of applicability 
at this temperature region as the change of the electric 
condensate value becomes too large. Indeed the mass 
obtained from the Stark effect at T > 1.09Tc becomes 
smaller than the smallest mass in Figs. [7] and [51 indi- 
cating the breakdown of the OPE in the Stark effect. 
Hence, we emphasize that any extrapolation of Fig. [TH 
to higher temperature is not appropriate. Apart from the 
marginal region, one sees that the downward mass shift 
smaller than the maximum given by QCD sum rules, as 
already found in Ref. [s^, leads to broadening just above 
Tc- These temperatures, 1.02Tc for rjc and 1.04Tc for 
J/V', are lower than the corresponding onset tempera- 
tures. One also sees that the continuum thresholds sud- 
denly decrease around Tc as in the case for the masses. 
It is quite intriguing to see that similar results are ob- 
tained in the potential model approaches, which utilizes 
the confinement force that can not be derived within the 
OPE formalism. Since the continuum thresholds change 
more rapidly, the resultant binding energy also drasti- 
cally decreases across Tc- At the marginal temperatures, 
i?bin is still around 100-200 MeV but the widths also be- 
come sizable due to thermal activation by gluons. We 
cannot draw conclusion on the dissociation of the char- 
monia from these results, since F > 100 MeV will have 
to be examined more carefully by incorporating higher 
dimensional operators and more realistic spectral func- 
tion. Below Tonsct, one finds -Ebin is still larger than F/2, 
indicating binding just above Tc- Furthermore, one does 
not see any broadening below T^. That is in line with 
our previous finding in Ref. [sll, where the effect of the 
continuum was not taken into account. Finally we would 
like to stress that all the spectral parameters show sud- 
den change across Tc, as shown in Fig.[I3 reflecting the 
abrupt change of the gluon condensates at this temper- 
ature and thus the QCD phase transition. Moreover, as 
discussed before, even if one of these parameters remains 
constant and retains its vacuum value, the QCD sum rule 
constraints force other quantities to exhibit such critical 
behaviors. 
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TABLE III: Spectral parameters of bb systems at T = 0. Ex- 
perimental masses are taken from Particle data book 59]. 

System ^ [GeV] m [GeV] m^xp [GeV] Mj [GeV^] /p [GeV^] 



T 

Xbo 

Xbi 



10.28 
10.34 
10.73 
10.34 



9.392 
9.447 
9.949 
10.09 



9.389 
9.460 
9.859 
9.893 



12.31 
13.68 
13.08 
14.13 



2.199 
2.034 

0.8 
0.492 
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FIG. 13: (color online). Same as Fig. |6] but for bb systems 



C. Results for bottomonia 

We also calculate the in-medium changes of the spec- 
tral property of the bottomonia using the same frame- 
work. Results for T = are summarized in Table IIIII 
As in the charmonium case, the sum rule works well for 
the bottomonium masses in the vacuum. In the case of 
bottomonium, the relative contribution of the dimension 
four operator to the OPE is much smaller than that of 
the charmonium because of the m7^ dependence in the 



TABLE IV: Onset temperatures of the width Tonsct for bot- 
tomonia 







T 


Xbo 


Xbi 


^onsct 


2.40 


2.56 


1.87 


1.50 
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[GeV] 

FIG. 14: (color online). Same as Fig. [71 but for rft- 
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FIG. 15: (color online). Same as Fig.0 but for T. 
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FIG. 16: (color online). Same as Fig. [T] but for Xbo- 



FIG. 17: (color online). Same as Fig. [71 but for Xbi- 



Wilson coefficient as seen in Eqs. ([T3| and Hence, 
its spectral property is much less affected by the change 
of the gluon condensates coming from the temperature 
effects. This fact allows us to go to much higher temper- 
atures than in the charmonium cases until the dimension 
four contributions become so large as to break the Borel 
stabilities. We display some of the OPE coefficients at 
T = and T = l.ST^ in Fig. One sees that the power 
correction terms are much smaller than the leading per- 
turbative correction term and might in fact be similar 
in magnitude to that of the next higher order radiative 
correction [53]. Therefore, while the separation scale in 
the present case is large enough for the OPE to provide a 
qualitatively reliable guide, further efforts are needed to 
obtain a quantitatively accurate estimate of the spectral 
property. ^From Fig. [T31 one also sees that the twist-2 
term dominates the temperature effect at 1.5Tc and also 
at temperatures where the spectral modification becomes 
sizable as will be seen below. Since the perturbative ef- 
fects are more dominant in the present case than in the 
charmonium detailed comparison with the re- 

summed perturbative approach j503 might be useful to 
understand the interplay between the perturbative and 
the non-perturbative effects at these temperatures. 

We show the constraints among the effective contin- 
uum threshold, mass and width for each of the bottomo- 
nium states in Figs. [TiHTTl The basic features are the 
same as in the charmonium cases, except now the sud- 
den change across Tc has disappeared. As the position 
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FIG. 18: (color online). Temperature dependence of the spec- 
tral parameters of T (closed symbols) and 776 (open symbols) 
extracted from QCD sum rules combined with the second or- 
der Stark effect. 



of the mass is far from 2mf, = 8.24 GeV, the unphysical 
behavior of the width seen at T > Tonsot in the charmo- 
nium cases is absent in the bottomonium cases. The on- 
set temperatures are summarized in Table IIVI One sees 
that the 5'-wave states have the narrow pole solution up 
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to T ~ 2ATc and P-waves do so up to 1.5Tc, suggesting 
survival of these states up to somewhat higher tempera- 
ture than those suggested by a potential model analysis 
[l3| . The maximum mass shifts are obtained near the 
onset temperatures and are found to be around 200 MeV 
for all the channels except for the axial vector. The mass 
shift of the P-wave states is found to be twice as large as 
that of the S'-wave states at a fixed temperature as in the 
charmonium states. Unlike the charmonium, one sees a 
significant difference in Tonsct between the S and A chan- 
nel. This is a manifestation of the dependence of Tonsct on 
the Borel window that was discussed in Sec. Ill CI In this 
case, the additional condition, the perturbative correc- 
tion less than 0.3 leads to a significant difference of M^j^ 
between S and A channel such that the Borel window of 
A channel becomes much narrower. If one set a common 
M^jjj for instance, Tonsot does not differ so much. Indeed 
the Borel stability in the Xb states is more marginal than 
other cases, as indicated in Ref. as the difficulty of 
establishing the "plateau" in the moment sum rule. For 
the A channel in the present case, even at T = 0, M^^^ 
is larger than Mg at ^/sq = oo. Although one can ob- 
tain the stability by relaxing the criterion for the M^^^, 
the result would be less reliable since the M^^^ is so cho- 
sen as to validate the perturbative expansion. After all, 
while we could obtain reasonable description of Xb states 
at r = and plausible in-medium changes of them, it 
has an intrisic ambiguity in the quantitative results. 

For the S'-wave states, we also extract the results from 
combining the constraints with the second order Stark 
effect which are expected to be more reliable in the bot- 
tomonium systems. Figure [18] shows the results for the 
mass, the effective continuum threshold, the binding en- 
ergy, and width of T and rji,. One sees that the changes as 
a function of the temperature are rather moderate; this 
reflects the smaller effects from the gluon condensates. 
Especially there is no significant brodening in both chan- 
nels up to T = 2.2Tc- One should note, however, that 
the second order Stark effect gives the larger mass shifts 
than that of maximum given by the QCD sum rule at 
T — 2ATc in V channel. This is similar to what hap- 
pened at T > 1.09Tc in the case of J/t/. (See Sec. iniB|. 
The heavier quark mass enables us to extend the OPE 
to higher temperature, but it seems to break down at 
T — 2ATc- Since the P channel exhibits larger spectral 
change than the V channel as in the charmonium cases, 
the maximum mass shift of rj^ given by the QCD sum 
rule is always smaller than that from the second order 
Stark effect. One sees small broadening at T > 2.2rc- 



IV. IMAGINARY TIME CORRELATORS 

The spectral parameters obtained from QCD sum rules 
at finite temperature have shown sizable modifications 
from the vacuum values. To confirm the findings, it is 
desirable to compare the results with the first principle 
lattice calculation. Unfortunately the direct evaluation 



of the spectral function of the heavy quarkonia through 
MEM has insufficient resolution to identify the spectral 
changes of order of 100 MeV. In this section, we will con- 
struct model spectral functions at finite temperatures as 
well as in the vacuum using the previously obtained QCD 
sum rule results. Then, we reconstruct the imaginary 
time correlators via the dispersion relation, discuss how 
the spectral modification affects the correlator, and then 
compare them with the lattice results which are more 
accurately calculated. 



A. Relation of spectral function with the 
imaginary time correlator 

The imaginary time correlator G{t, T) is related to the 
spectral function via the dispersion relation 



G(r,T) 



dujK{uj,T;T)p{uj,T) (21) 



where the integration kernel K{uj,t;T) is 



K{uj,T-T) 



cosh[u;(T- l/2r)] 
sinh(w/2r) ' 



(22) 



of which the zero temperature limit is e~"'^. To see the 
temperature effect on the spectral function, one usually 
computes the ratio of this correlator to the reconstructed 
one G{t,T)/Gi-cc{t,T) with Gicc defined as 



/>oo 

Grcc(T, T)= / (hjpiLJ.T 

Jo 



0)K{lj,t;T) (23) 



which has temperature dependence coming only from the 
kernel. 

We construct a model spectral function to be put into 
Eq. ((2T|) from the phenomenological side p^ - pT)) 



pP^(c.) 



ImIIP°^"(a;2) + Imn'=°"'(w2)l (24) 



with Gj = 1 for P and S channels and 3 for V 
and A channels. The subscript "pc" denotes the 
"pole-|-continuum" . We relate them to the spatial com- 
ponents of the spectral function for V and A channel, 
in order to compare them with lattice calculation. In A 
channel, although lattice calculation uses an axial vector 
current of = hj^j^h while we use the conserved part 
jAi ^ jfj^ by muhiplying ijap = [qaqp/q^ - ffa/s), the 
above expression still holds. 

Putting each part of the model spectral function into 
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Eq. ([2T|) . one obtains the following formulae; 

' Cjmfo cosh[m(T - 1/2T)] 



GP°"=(r,T) 



27r 
CjfV 



sinli(m/2r) 



cos1i[lj(t- 1/2T )] 
sinli(a;/2r) 



r = 



r ^ 

(25) 



TABLE V: Numerical constants in zero mode spectral func- 
tion of various channels. 
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S yl(p") 
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C2 
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2 


-2 -4 


-6 



G™"*(t.T) 



cosh[cj(T - 1/2T)] channels are given by 

p'--\Lo)^N,uj6{u:){cih-C2l2)- 



sinh(w/2r) 
(26) 



(29) 



As for the peak strength parameter / and /g, One can 
obtain it by using the dispersion relation (|15p after de- 
termining the other three parameters as 



/o = e™'/*^" [X(M2) - 7W™"*(m2)], 



r 



ds f 



(27) 
(28) 



Af2 



We adopt the value of the Borel mass at — 
where the property of the pole part was determined. This 
choice, however, can result in an unphysical behavior of 
the reconstructed correlator ratio G/Grec near r ~ due 
to the sensitivity to the high energy part of the specral 
function. The continuum part in the P and the S chan- 
nel has a singular behavior in the high energy limit (ioj 
such that a slight deviation in as can lead to a sizable 
difference in the imaginary time correlator near r = 0. 
While this does not matter at r > 0.1 fm and in the 
QCD sum rule analysis due to the large suppression of 
the high energy part by the Borel transformation, we use 
the same the value of as at finite temperarature as that 
of T = by fixing so that G/G^cc -J> 1 as r at 
any temperature. 

Explicit temperature, not divided by Tc, need to be 
specified in the kernel. While our gluon condensates have 
been taken from the lattice calculation with Tc — 264 
MeV (42|, we normalize the temperature dependence in 
the imaginary time correlator calculation to Tc — 295 
MeV, which corresponds to the normalization used in the 
lattice calculation that we will be comparing our results 
to 



22| 



It has been emphasized that a peak of the spectral 
function at w = gives a constant contribution to the 
imaginary time correlator [45| . Although we have ignored 
this contribution in the QCD sum rules, as explained 
above, this term is necessary for proper comparison of 
G(t, T). Here, we adopt the expression calculated for free 
heavy quarks which is proportional to ujS{lj). In this case, 
the spectral functions have been calculated and given in 
[ssj.^ The zero mode (scattering) parts for V, P, S and A 



There is a misprint in Ref. ISSH pointed out in Ref. 



The numerical constants ci and C2 are summarized in 
Table |V] and 



h = 



-2 



d^k duk 
(27r)3 doJk 

d^k duk 
(27r)3 dujk L^l 



(30) 



where — (e"^*"/^ -|- 1) ^ and ujk = y k + mf. Putting 

these expressions into Eq. (PT|) . finally one obtains the 
constant contribution to the imaginary time correlator 



G^--\T,T) = NcTicih-C2l2). 



(31) 



Hereafter, we adopt the three component model 
G(r, T) = GP°i'= + G™"* -I- G"'^''* with spectral parame- 
ters taken from the results of the QCD sum rule as our 
model imaginary time correlator which we compare with 
the lattice QCD result shown in Ref. 22]. 

Before proceeding, it is useful to see how the typical 
spectral changes seen in the QCD sum rule affect the ratio 
of the imaginary time correlator G/Groc- In Fig. [111 we 
plot several examples in which only one of the four spec- 
tral parameters is changed. Similar to what was shown in 
subsection IIII A[ one can understand the qualitative be- 
havior from the dispersion relation for the imaginary time 
correlator (|21l) due to the positivity of the kernel. When 
the modification of the spectral parameters increases the 
spectral sum, G/Grec also increases. One sees that a 
rather small modifications of the continuum threshold 
and the mass, of order of 100 MeV, lead to more than 
10% change in G/G^cc, whereas less than 3% changes 
have been observed in the lattice QCD calculation of the 
P channel 22]. On the other hand, the change of width 
does not affect G/Groc, as seen in the lower- right panel. 
This is so because the increase of width implies reduc- 



tion of the peak height p{uj 



^) fm/T. If one 



artificially tries to preserve the height by increasing / as 
well as r, even small increase of the width makes G/Grcc 
increase very quickly, as shown in the dotted line in the 
figure. 
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FIG. 19: (color online). G/Gmc in several cases of modifi- 
cation of single spectral property. Upper left : continuum 
threshold, upper right : mass, lower left : width, and lower 
right : peak height. The zero temperature spectral function 
in Groc is taken from rjc case thus constant contribution is 
absent. For the lower panel, we utilized T — spectral pa- 
rameters obtained by putting F = 2 MeV for comparison. 

B. Comparison with lattice QCD results 

1. Charmomum 

We start with the charmonium in the P channel (rjc) 
to which no zero mode contributes. Unfortunately, the 
available lattice data are for T = 0.87Tc, l.OTTc, 1.09Tc 
and so on while our sum rule results of interest are just 
around Tc. At O.STTc, spectral change is almost negligible 
because of the tiny change of the gluon condensates. At 
the next lowest temperature, 1.07Tc, it is already above 
the onset temperature therefore our result for the char- 
monium is quantitatively obscure. Between the lowest 
and the next lowest temperature, however, G/Groc on 
the lattice seems unchanged because it is almost unity 
at both temperatures. Sizable deviation from unity has 
been seen above l.STc at which the spectral functions ex- 
tracted by MEM also show notable modifications. Hence, 
for references, we calculate G/G^ec at not only tempera- 
tures where the lattice QCD results are available but also 
at r = Tc and at T = IMT^ where our QCD sum rule 
works well and results of the lattice QCD are unambigu- 
ously anticipated. 

Figure [201 displays G/Groc for the P channel charmo- 
nium current. We plot the cases for several sets of the 
spectral parameters, summarized in Appendix [Cj At 
T = Tc, G/Grec of rjc can be both below and above unity 
depending on the parameter set. The parameter set of 
=3.4 GeV which seems most plausible according to 
Fig. [T2] gives the closest result to unity. We anticipate 



FIG. 20: (color online). Imaginary time correlators for rjc. 
Each of panels shows different temperature cases, T — Tc, 
IMTc, 1.07Tc, and 1.09Tc. For T = imTc and T = 1.09Tc, 
we plot results from lattice QCD shown in Ref. [2^ ] with 
crosses. 
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FIG. 21: (color online). Integration kernels of the dispersion 
relations for the correlation functions. 



the lattice QCD gives almost unity as in the higher tem- 
perature cases. Therefore the most plausible case gives 
the most consistent result with the lattice result. Quali- 
tatively, our results indicate that the combination of the 
small decreases of both effective threshold and mass can 
give a consistent result. One sees, however, that all of 
the parameter sets lead to G/Grcc > 1 above T > 1.04rc 
in spite of the quite different spectral parameters. We 
would like to stress that T = 1.04rc is the onset tem- 
perature of rjc] for the parameter set with the smallest 
y/so shown in the figure, we need to introduce F = 18 
MeV to stabilize the Borel curve. No Borel window is 
open for smaller y/sQ. In the two higher temperature 
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cases, T = 1.07Tc and 1.09Tc, lattice data are shown to- 
gether. Note that the lattice setup is different between 
these temperatures. The latter is obtained with finer 
lattice spacing. As shown in the figure, the lattice cor- 
relators do not exhibit sizable changes from unity, that 
would be interpreted as no spectral modification at these 
temperatures. 

These disagreements with the lattice QCD results at 
T > 1.04rc do not immediately mean serious flaws in 
our approach. First, the applicability of our present ap- 
proach to these temperatures are marginal and thus our 
results will be quantitatively improved by including cor- 
rection mentioned above. Second, our model spectral 
function might be too simple to make such a compar- 
ison. This simplification does not matter in the QCD 
sum rule approach due to the strong suppression of the 
high energy part in the Borel transformed dispersion re- 
lation, but might cause this defect in the correlator be- 
cause of its sensitivity to the continuum. Figure [3T] dis- 
plays the kernel in the dispersion relation (1221) with typ- 
ical values of parameters together with the kernel in the 
Borel transformation (jlSp where we rewrite the formula 
so that integration variable is a; = ^/s and integrand has 
a form of K{uj, ]\'P)p{ijj). One sees that the Borel trans- 
formation suppresses the high energy part of the spectral 
function much more strongly than the heat kernel (1^^ . 
Note that temperature effects on the heat kernel are small 
around T^, it almost behaves like e"""^. This fact means 
that while the detailed structure of the continuum and 
the excited states do not affect the property of the pole 
part in the QCD sum rule approach, they do contribute 
non-trivially to the imaginary correlator. For example, 
compared to the realistic situation, we have neglected 
the contributions from the excited states (77c (25") in P 
channel) in the model spectral function. While such ap- 
proximation does not affect the lowest pole in the QCD 
sum rule analysis, it will result in the smaller continuum 
threshold than a realistic value to compensate for the 
missing state. This discrepancy between effectively small 
continuum threshold and existence of excited states will 
result in different imaginary time correlators. If the 25* 
state melts just above (or below) Tc, our parametrization 
is in fact better above Tj,. Hence, if one takes into ac- 
count the excited state contribution explicitly, Grec (t, T) 
will increase while G(t, T) remain unchanged, thus the 
ratio will now become close to unity. Indeed the effect of 
the excited state has been examined in the context of a 
temperature dependent potential model in Ref. [ll| ■ The 
authors found that the r/c(25) state reduces G/Gjcc 10- 
20 % because the 25* state melts above Tc thus increases 
Gicc- One sees that the observed reduction is also rele- 
vant for our case. Although the result of Ref. [11] shows 
the 2S states melts at Tc thus our agreement at this tem- 
perature might be changed to other set of parameters, the 
spectral property of the excited states at finite tempera- 
ture is not clear yet. If the 25' state survives at Tc and 
our analysis catches essential points for the correlator, 
the difference between T = Tc and T > Tc indicates that 
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FIG. 22: (color online). G/Grcc for the V channel charmo- 
nium current. 

the 25" state dissolves at < T < l.MTc. 

In Fig. [221 we display the results of G/G^cc for the 
vector channel (J/ip)- One sees a deviation from unity 
starting at small r region. This behavior is caused by 
the zero mode contribution [13, 45, sj]. At T = 1.07rc, 
our G/Groc with ^fso < 3.2 GeV agrees the lattice results 
at T < 0.15 fm in which the zero mode contribution is 
negligible. This means that G/Gicc — 1 is achieved by 
a combination of the various spectral changes as seen in 
the case of r]c at T = Tc- The zero mode contribution, 
however, overwhelms other changes at T = 1.07rc and 
T — 1.09Tc, as clearly seen from the figure. There are 
following possibilities; 

1. High energy part of the model spectral function. 

As in the case of i]c, we have neglected ip'{2S) con- 
tribution to the T — spectral function. Includ- 
ing this leads to larger Groc thus reducing G/Groc- 
There may be also the possibilities that we under- 
estimated the pole modification by truncation of 
the OPE and other approximations. As for d = 6 
contribution in the OPE, however, it is expected to 
reduce the spectral modification |57| . 

2. Free charm quark approximation in the zero mode 
contribution. 

Indeed this might be a flaw because in more re- 
alistic situation quarks are interacting such that 
zero mode spectral function is smeared 5S5- Of 
course, since we are looking at the integrated value 
of the spectral function, this smearing itself might 
not change the value so much. Nevertheless, there 
are further ambiguities in the zero mode calcula- 
tion such as the value of the charm quark mass; 
within the quasiparticle picture, the thermal effect 
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FIG. 23: (color online). Effect of 25 state in the imaginary 
time correlator. Upper panel : ratio of the correlator G with- 
out -i/)' to that with ip' a,t T — 0. Lower panel : G/Grec 
at T = l.OTTc with tp' contribution being included in Groc. 
Thin lines are the same ones shown in Fig. [22] while thick lines 
denote the cases in which tp' contribution is included in Grcc . 

will effectively increase the quark mass such that 
the zero mode contribution is reduced according to 
the thermal distribution. As we shall see below, 
there is clearly something that cannot be under- 
stood with the free charm description in the zero 
mode. To avoid these ambiguities, subtraction of 
the zero mode contribution by taking derivative the 
imaginary time correlator and then looking at the 
ratio [isl |45| will provide useful information. 

In the vector channel, the first possibility can be ex- 
plicitly checked by including ip' contribution to the con- 
tinuum part of the model spectral function at T — 0. 
Namely, 

lmU^'is) = fS{s-ml,), (32) 

with /' = 0.276 GeV^ being obtained from the lep- 
tonic decay width given by the Particle Data Group 
[59| , is added to the phenomenological side ([TC)) . With 
this implementation, we found that the resultant J/ip 
mass changes only 0.3% (3.05 GeV) while the contin- 
uum threshold increases from y^=3.54 GeV to 3.93 
GeV. On the other hand, we also found that incorpo- 
rating the ip' to the spectral function in the dispersion 



relation of the imaginary time correlator ([2T|) leads to a 
sizable change. This fact exactly demonstrates our ex- 
pectation discussed above; while the spectral property 
deduced from the QCD sum rule is independent of de- 
tailed structure of the higher energy part of the model 
spectral function as long as the pole dominance is well 
satisfied, the imaginary time correlator receives sizable 
change. Note that the physical meaning of the thresh- 
old parameter becomes different if one takes the excited 
state into account. When one includes it, now can 
be regarded as physically more relevant threshold while 
it represents an effective one controlling the contribution 
from the high energy part other than the lowest pole in 
the dispersion relations. For example, y/so in Fig. [12] 
means a merely effective threshold parameter since we 
have neglected the excited state contribution through- 
out the calculation. If ip' melts above Tc, can be 
now regarded as more physical one. In this case the true 
threshold might vary more rapidly from 3.93 GeV to a 
lower value in the vicinity of Tc- 

Figure [53] shows the effect of ip' on the imaginary 
time correlator. In the upper panel, we compare two 
G{t,T — 0), with and without -0' contribution. One sees 
3% reduction of the ratio, which means that inclusion of 
the ip' gives a enhancement large enough to affect G/Grcc 
comparison. The resultant G/Grec are shown in the thick 
lines in the lower panel, where one sees the apparent re- 
duction of G/Groc when including ip' . Nevertheless, r 
dependence is still governed by the zero mode contribu- 
tion, on which we will give further consideration. 

As for the second possibility, one way to check the 
consistency is to look at other channels. We display the 
imaginary time correlators of the scalar and the axial- 
vector channels in Figs. [M] and [25] respectively. Simi- 
larly the spectral parameters are summarized in Table 
IVIIII and IIXI These two channels show quite similar be- 
havior so that the following discussion can be applied 
for both cases. First, no sizable difference among vari- 
ous parameter sets is seen as indicated by the complete 
overlaps of the lines. One sees the clear effect of the 
zero mode contribution and its agreement with the lat- 
tice results at T = 1.07Tc contrary to the V channel case. 
This might be partly attributed to the absence of the 25' 
state contribution below the continuum threshold in the 
S and the A channels; i.e., the single pole plus continuum 
ansatz at T = is a better approximation in these chan- 
nels than in the V and the P channels. At T = 1.09Tc, 
however, this agreement is lost though qualitatively the 
lattice results indicate the dominance of the zero mode. 
One sees the value of G/Gioc is smaller at T = 1.09rc 
than at T = l.OTTc in the lattice results. If the spectral 
modification of the pole and the continuum part does 
not differ so much between these temperatures, this re- 
sult seems to indicate the smaller zero mode contribu- 
tion at higher temperature. This cannot be understood 
within the free charm quark approximation in which zero 
mode contribution increases as temperature increases if 
the charm quark mass is constant. Therefore, although 
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FIG. 24: (color online). Same as Fig. [20l but for 5 channel 
current. 
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FIG. 25: (color online). Same as Fig. 1201 but for A channel 
current. 



FIG. 26: (color online). Ratio of the derivative of the imag- 
inary time correlator G'/G^ec for the V channel. Shown pa- 
rameter sets are the same as in Fig. 1221 
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FIG. 27: (color online). Same as Fig. 1261 but for S channel. 



our results show agreement with the lattice results in the 
T range where zero mode contribution is relatively small, 
we cannot draw definite conclusion on the quantitative 
correctness of the zero mode from the results. One 
way to avoid the difficulty of the zero mode is to evalu- 
ate the derivative of the imaginary time correlator with 
respect to r [i^ . Figures USHM] show the ratio of the 
derivative of the imaginary time correlator, G' /G[^^, in 
which the zero mode contribution is absent. One sees 
different tendency from G/Grec such that G'/G[^^ ~ 1 
within 3% for certain sets of the spectral parameters at 
T < Tonset- This strongly supports robustness of our 
results at these temperatures since the lattice computa- 



tions, although not available at these temperature, are 
expected to be unity also. We would like to note that 
the parameter set close to the ones constrained from the 
Stark effect gives G' /G'^^^ closest to unity in the case of 
J/ip- For illustrative purpose, we display the spectral 
density of J/ip constructed from Eq. in Fig. with 
the same parameter sets as those in Fig. [T^l 

As demonstrated in Ref. [11], effect of the threshold 
enhancement [l^ . [isj seems important for understand- 
ing the relation between the lattice measurement which 
gives G/Gicc — 1 and model calculations such as poten- 
tial models. In rjc and J/ip, it explains successfully the 
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FIG. 28: (color online). Same as Fig. 1261 but for A channel. 
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FIG. 30: (color online). G/dcc for bottomonium currents. 
Lattice results taken from Ref. [131 . Left and right columns 
stand for T — l.lSTc and T = l.SiTc cases while each rows 
denote P, V , and S channels from top to bottom, respectively. 



inance will be badly violated, this behavior of y may 
indicate a connection to the higher temperature regime 
which has been investigated only through the potential 
models so far. 



2. Bottomonium 



FIG. 29: Spectral density of J/^l^ at the best fit to the lattice 
correlator. 

lattice data at T = 1.2rc. Note that this is beyond our 
regime in the present work; we cannot apply our method 
at this temperature since we could not have the reliable 
Borel stability beyond Tonsct. We can still see some in- 
dication of such an effect in the present data. One notes 
that the Borel curve becomes flatter and flatter, accord- 
ing to ^/x^, as decreases above T = Tonsct. (See 
Tables IVTl and IVll|) . It has minimum at the smallest 
at which the Borel window is about to close. This fact 
means that if one relaxes the criterion for the pole dom- 
inance, one still obtains the stability for smaller y/sQ. 
Though any extrapolation to higher temperature can- 
not be reliable due to the missing effects as discussed, 
one might expect the threshold to be close to the mass 
thus the spectral function exhibits the threshold enhance- 
ment. Whereas direct confirmation will not be possible 
with the present method by construction since pole dom- 



We compute G/Groc for the bottomonium currents 
in the same manner. Figure [30] displays the results at 
T = l.lSTc and T = IMTc of which lattice results are 
available with the finest lattice spacing in Ref. [12]. As 
in the charmonium cases, the lattice results show lit- 
tle deviation from unity. In the 5- waves, while our re- 
sults show variations among the parameter sets there is 
a certain range of the effective threshold parameter of 
which parameter set gives G/Grcc — 1. The best agree- 
ment seems to lie between the data set with vacuum 
(green, long-dashed) and that with Stark effect results 
(blue, short-dashed). The difference of these two cases is 
tiny at T = l.lSTc but sizable at 1.54Tc, indicating the 
possibility of discriminating the spectral change from the 
imaginary time correlator. As before, one has to consider 
excited states to give a definite conclusion. The interpre- 
tation of the behavior of G/Groc depends on whether such 
states below threshold survive or not at these tempera- 
tures. From these agreements, a possible interpretation is 
that excited state are still surviving even at T = 1.5ATc 
since a potential model calculation shows larger reduc- 
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FIG. 31: (color online). G'/Gj-oc for bottomonium currents. 
Left and right columns stand for T = 1.15Tc and T = 1.54rc 
cases while each rows denote V , S, and A channels from top 
to bottom, respectively. 



tion of G/Gicc by including 25* and 3S states in rji, than 
the variation seen in the figure 11]. The result of Xbo 
state is again dominated by the zero mode contribution 
as was the case for the Xco- One sees that even lattice 
results show a difference between that at T = 1.15Tc 
with /3 = 6.1 and that at T = l.lGTc with /3 = 6.3, 
indicating the difficulty in quantifying the calculations. 
Therefore we compute G' /G[^^ as was done in the char- 
monium cases. Figure [31] shows the results for G' /G[^^ of 
V, S, and A channels. As before, one sees that G' /G[^^ 
is sensitive to the variation of the spectral parameters 
and there exist certain ranges which give G' /G[.^^ ~ 1. 
Precise determination of this quantity will be useful for 
constraining the spectral changes. One sees that drop- 
ping of both the mass and the continuum threshold is 
consistent with G' /G[.^^ ~ 1 and that the result from the 
second order Stark effect again fits well in the case of T. 



V. SUMMARY AND OUTLOOK 

We have analyzed the spectral changes of heavy 
quarkonia in the hot environment in a systematic way 
based on the QCD sum rule with Borel transformation 
technique. We have taken into account possible changes 
of the continuum spectrum by a temperature dependent 
effective continuum threshold. Although the tempera- 
ture dependence of the OPE side allows various 
combination of the changes of the spectral parameters. 



we have given the constraints among them by an opti- 
mization procedure which has been widely used in QCD 
sum rule applications. We found that instability of the 
Borel curve in the F = limit caused by the change 
of the gluon condensates can be cured by introducing a 
width although neglected contributions from higher di- 
mensional operators might be non-negligible at the small 
Borel mass. The results, FigslTHTUl and [T1HT71 show the 
behaviors of the spectral parameters with respect to the 
change of temperature through the gluon condensates. 
As already argued in previous literatures [13, [H, [39| . 
charmonia exhibit the critical behavior in at least one of 
the spectral parameters. Note that although the effective 
continuum threshold shares the effect from the change of 
the gluon condensates and thus reduces the change of 
other parameters, it is basically linked with the mass as 
seen in the Figs. [7l-[T0l and [MHTTI as a result of stabiliz- 
ing the Borel curve. While we do not taken into account 
the excited state explicitly, property of the lowest pole 
is not affected by this simplification since it is imposed 
on the effective threshold parameter. When one of the 
spectral parameters remains unchanged, rapid change of 
other parameters is inevitable. We found there is an on- 
set temperature for each channel at which broadening 
must occur. Although the value of this temperature is af- 
fected by the assumption on the Borel window and does 
not necessarily mean there is no significant broadening 
below, the combined analysis with the QCD second or- 
der Stark effect shows that r/c and J/ V' do not likely have 
significant width below Tc and all the spectral parame- 
ters of J/ip and r]c change abruptly in the vicinity of T^. 
The same analysis procedure for the bottomonia shows 
the little change around Tc because of the much heavier 
quark mass, but eventually shows sizable changes with 
increasing temperature. 

Such changes should be obtained from the lattice QCD 
also, while spectral analyses based on MEM do not have 
sufficient resolution. Therefore, we have computed the 
imaginary time correlator which is the basis of the MEM 
analyses and its derivative with respect to the Euclidean 
time r by putting the phenomenological side of the QCD 
sum rule analyses as a simple model of the spectral func- 
tion. Then, we take the ratio G/Groc in order to see 
the temperature effect on the spectral function, as done 
in the lattice analyses and potential model calculations. 
We have demonstrated that the results obtained in the 
lattice calculation, namely G/Gioc — 1, do not always 
mean absence of the spectral changes, but can mean a 
mixture of some sets of the change of the spectral proper- 
ties. Similar observation has been done in some potential 
model analyses [H, [l3| but they are at higher tempera- 
tures and mainly focused on the threshold enhancement. 
We showed that the rather small modification compared 
to the bound state masses gives the sizable change in 
G/Grcc and the constrained parameter sets by the QCD 
sum rule can lead to G/G^ec — 1- We have also pointed 
out that while similarity in the dispersion relation exists 
between the Borel-transformed current correlation func- 
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tion and the imaginary time correlator ([2T|) , the for- 
mer is much more dominated by the pole contribution. 
As a result, the QCD sum rule analysis is not much af- 
fected by taking into account the known excited states 
explicitly, while the imaginary time correlator shows a 
small but significant change over the uncertainties of the 
lattice QCD results. Due to this property and the poorly 
known spectral function near 2mc threshold region in the 
case of finite F, more work is needed to give a precise 
quantitative determination of the spectral parameters. 

Before closing, we comment on implication for the full 
QCD case by repeating our argument in Refs. (27l[28l.l39| 
since the present work is based on quenched QCD. If one 
takes into account light dynamical quarks in the OPE, 
they appear as light quark condensate contributions and 
change of the temperature dependence of the gluon con- 
densates in the OPE. The former can be safely neglected 
since it is at order aKq^). The latter has been shown 
to lead smoother temperature dependence of Go near Tc 
[13] reflecting the crossover nature of the transition. One 
also expects the similar change in G2- Since the OPE 
side depends on temperature through these condensates, 
this will result in smoother spectral changes than those 
shown in the present work. While the spectral changes 
become smoother, the actual magnitude of changes at Tc 
might not differ so much since the reduction of the scalar 
condensate at T = Tc is almost the same as that for the 
quenched case [27]. This also implies that the broaden- 
ing below Tc could be negligible. On the other hand, the 
moderate decrease of the condensates above Tc may lead 
to higher onset temperatures in the full QCD case. 

The results shown above strongly indicate that our 
main results, the mass shifts and width broadening in- 
duced by the QCD phase transition, might be also real- 
ized in lattice QCD simulation. The agreement between 
the sum rule constrants and the lattice correlator ratio is 
obtained on the basis of the effect of the known excited 
state on the lattice correlator and the zero mode con- 
tribution of the free heavy quarks. The best set of the 
parameters depends on such external assumptions. Fur- 
ther assessments of these quantities as well as continuum 



spectrum at finite temperature [5(| will be required to 
check the consistency more quantitatively. 
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Appendix A: Borel transformed Wilson coefficients 

In this appendix, we list the Wilson coefficients seen in 
Eq. (fT^ originally obtained by Bertlmann in Ref. pl| . 
Hereafter 
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Appendix B: Spectral function for the continuum part 



We summarize the continuum part of the phenomenological spectral function Imn'^'P'^''*(s) which are taken from 
the perturbative QCD calculation up to 0{as) shown in Ref. j3^]. Here, 
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Appendix C: List of obtained spectral parameters 

We summarize the results of the spectral parameters of 
the model spectral function at various temperatures for 
which we compared them with the lattice QCD results of 
the imaginary time correlators in Sec. IIVI In the case of 
r = 0, the minimum ^ will provide a good estimation 
of the best Borel curve since M^jj^ is fixed. When F > 0, 
however, this will no longer hold to choose the best one 



among curves with different ^/sq. Therefore, we used this 
criterion only to choose the best F value among a fixed 
case for each temperature. In Tables IVlllIXl we list 
sets of resultant parameters obtained in this way (not all 
cases: we computed more cases of ^sq but we include the 
case of the smallest ^Jsq below which no Borel window is 
available). We list a/x^ instead of since it is related 
to a crude estimate of a systematic error in the mass as 
discussed in Sec. Ill CI 
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TABLE X: Summary of spectral parameters for T obtained from QCD sum rules. 



T/T, 


[GeV] 


m [GeV] 


r [MeV] 


/ or /o [GeV^j 


[GeV^] 


[MeV] 


1.15 


10.04 
10.16 
10.28 
10.34 
10.44 


9.369 
9.401 
9.428 
9.440 
9.480 







40 


1.670 
1.845 
1.959 
2.023 
0.694 


16.52 
14.86 
13.71 
13.26 
13.31 


44.3 
31.7 
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24.9 
29.3 


1.54 


10.0 
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10.2 
10.34 
10.5 
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9.377 
9.401 
9.441 
9.516 
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1.613 

1.739 
1.857 
0.647 
0.738 


16.20 

14.74 
13.68 
13.03 
13.46 


38.3 
27.9 
21.6 
23.0 
30.3 




TABLE XI: Summary of spectral parameters for rji, obtained from QCE 


1 sum rules. 






^ [GeV] 


m [GeV] 


r [MeV] 


/ or /o [GeV^] 


[GeV^] 


V3? [MeV] 


1.15 


9.98 
10.1 
10.22 
10.28 
10.4 


9.315 
9.347 
9.373 
9.385 
9.430 
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1.769 

1.941 
2.094 
2.163 
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12.33 
11.91 
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48.6 
35.1 
28.2 
28.3 
34.7 
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9.94 
10.04 
10.14 
10.28 
10.4 


9.295 
9.322 
9.345 
9.387 
9.441 






28 
86 


1.705 
1.848 
1.977 
0.694 
0.768 


15.18 
13.23 
12.17 
11.69 
11.95 


39.5 
28.9 
22.9 
25.7 
31.8 




TABLE XII: Summary of spectral parameters 


for Xbo obtained from QCD sum rules. 




T/T, 


^ [GeV] 


m [GeV] 


r [MeV] 


/ or /o [GcV'l 


M'^ [GeV'l 


[MeV] 
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10.73 
i().8 
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9.98 1 
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10 
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I2.8:i 


22.5 

13.4 
16.5 
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9.839 
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iO.OOT 



16 

84 

110 


0.542 
0.207 
0.260 

I). 28:! 


12.36 
11.89 
12.58 

12.8:! 


9.13 
7.17 
14.5 

18.1 


TABLE XIII: Summary of spectral parameters for Xbi obtained from QCD sum rules. 


T/T, 


[GeV] 


m [GeV] 


r [MeV] 


/ or /o [GeV^j 


[GeV^] 
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3.7 
5.4 
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